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ABSTRACT 

A three-dimensional direct numerical simulation (3D DNS) is performed to describe the turbulent flow in 
an enclosed rotor-stator cavity characterized by a large aspect ratio G = {b — a)/h = 18.32 and a small 
radius ratio a/b = 0.15 (a and b the inner and outer radii of the rotating disk and h the interdisk spacing). 
Recent comparisons with velocity measurements [1 ] have shown that, for the rotational Reynolds number 
Re = £lb 2 /v = 95000 fO the rate of rotation of the rotating disk and v the kinematic viscosity of water) 
under consideration, the stator boundary layer is 3D turbulent and the rotor one is still laminar. Budgets for 
the turbulence kinetic energy are here presented and highlight some characteristic features of the effect of 
rotation on turbulence. A quadrant analysis of conditionally averaged velocities is also performed to identify 
the contributions of different events (ejections and sweeps) on the Reynolds shear stress. 



Introduction 



Besides its primary concern to many industrial 
applications, the rotor-stator problem has proved 
a fruitful means of studying turbulence in con- 
fined rotating flows. This specific configuration 
is indeed among the simplest ones where rota- 
tion brings significant modifications to the tur- 
bulent field. The present paper is devoted to the 
study of the turbulent flow in an enclosed high- 
speed rotor-stator system of large aspect ratio. 
The flow is of Batchelor type and belongs to the 
regime denoted IV by Daily and Nece [2] . These 
authors provided an estimated value for the local 
Reynolds number at which turbulence originates 
with separated boundary layers, Re r = Vtr 2 jv = 
1.5 x 10 5 (r the radial location) for G < 25. 
However, experiments have revealed that transi- 
tion to turbulence can appear at a lower value 
of the Reynolds number within the stator bound- 



ary layer (the Bodewadt layer), even though the 
flow remains laminar in the rotor boundary layer 
(the Ekman layer). Recently, the 3D computed 
solution presented here has been compared to ve- 
locity measurements performed at IRPHE [1]. It 
has been shown that, for the rotational Reynolds 
number Re = 95000, the Bodewadt layer is tur- 
bulent and the rotor one is still laminar. The pur- 
pose of this work is to provide detailed data of 
the turbulent boundary layer along the stator side 
in a large enclosed system (G = 18.32). 



Numerical approach 



The numerical approach is based on a pseudo- 
spectral technique using Chebyshev polynomials 
in the radial and axial directions with Fourier se- 
ries in the azimuthal direction associated with a 
semi-implicit second order time scheme. An ef- 



ficient projection method is introduced to solve 
the pressure- velocity coupling. A grid resolution 
composed of N x M x K = 300 x 80 x 100 
respectively in radial, axial and azimuthal direc- 
tions has been used, with a dimensionless time 
step St = 2.75 x 10~ 3 . Numerical computations 
have been carried out on NEC SX-5 (IDRIS, Or- 
say, France). 



Results and discussion 



The detailed description of the mean field and of 
the axial variations of the Reynolds stress tensor 
is given in [1]. Nevertheless, we recall the main 
results. A good agreement has been obtained be- 
tween the 3D solution and the experimental data, 
whereas the axisymetric solution leads to a steady 
laminar flow. The flow is of Batchelor type: the 
two boundary layers developed on each disk are 
separated by a central rotating core characterized 
by a quasi zero radial velocity and a constant 
tangential velocity. The Bodewadt layer along 
the stator is centripetal, three-dimensional as the 
Townsend structural parameter is lower than the 
limit value 0.15, and turbulent with turbulence in- 
tensities increasing from the axis to the periphery 
of the cavity. On the contrary, the Ekman layer 
on the rotor is centrifugal and laminar. The tur- 
bulence is observed by the formation of turbulent 
spots developing along spiral arms towards the 
periphery of the cavity, as seen in figure 1 from 
the instantaneous iso-values of the turbulence ki- 
netic energy within the stator boundary layer. 



Turbulence kinetic energy budgets 



The balance equation for the turbulence kinetic 
energy k can be written: 




Fig. 1. Iso-contours of the instantaneous turbulence 
kinetic energy within the stator boundary layer. 

with the advection term A, the production term 
P, the turbulent diffusion D T , the viscous diffu- 
sion D u , the velocity-pressure-gradient correla- 
tion II and the dissipation term e. 
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Fig. 2. Budgets for the turbulence kinetic energy 
k/(nb) 2 at: (a) r/b = 0.56, (b) r/b = 0.8. 



A = P + D 1 +D U + n-e 



(1) 



Figures 2a and 2b show the axial variations along 
z* = z/h of the different terms involved in the 



transport equation (1) at two radial locations. It 
is clearly seen that all these terms vanish towards 
the rotor side (z* = 0), confirming the lami- 
nar nature of this zone up to the stator bound- 
ary layer (z* = 1). At the stator wall, the vis- 
cous diffusion balances the dissipation as well 
known in 3D turbulent boundary layer. Within 
the Bodewadt layer, even though some interac- 
tion between the different terms involved is ob- 
served, the major contributions come from the 
production, the dissipation and the viscous dif- 
fusion terms. The production is balanced by the 
dissipation and the viscous diffusion, the level of 
which increases in association with the thicken- 
ing of the boundary layer towards the periphery. 
The production increases with increasing radius 
as already observed with the levels of the nor- 
mal Reynolds stresses [1]. The maximum of the 
production term is obtained at the wall coordi- 
nate z + = zv T /v = 12 (v T = {{vdVe/dzf + 
(vdVr/dz) 2 ) 1 / 4 the total friction velocity and z 
the axial coordinate) for r/b — 0.56 and at z + = 
12.5 for r/b — 0.8, which confirms the approx- 
imately self- similar behavior of the Bodewadt 
layer. The levels of the viscous diffusion increase 
when moving towards the outer casing, where the 
highest turbulence intensities prevail. It indicates 
that viscous effects still play an important role in 
the turbulence towards these regions, which does 
not allow for a distinct delineation of the viscous 
sublayer. This indicates also the weak nature of 
the turbulence obtained at this Reynolds number. 

Conditional-averaged quadrant analysis 

To gain a better insight on the near-wall struc- 
ture of the turbulent boundary layer along the 
stator side, a conditional- averaged quadrant anal- 
ysis is performed to identify the contributions 
of intense intermittent bursting events (ejections 
and sweeps) on the Reynolds shear stress pro- 
ducing vortical structures. It corresponds to four 
subdivisions of the fluctuation field according to 
the combination of the tangential v' e and axial 
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Fig. 3. Variation with — Ar + of the conditionally 
averaged Reynolds shear stress at z + = 17 in the 
vicinity of (a) a strong ejection < v' g v' z \ strong 
ejection > / < v' e v' z > and (b) a strong sweep 
< v' e v' z I strong sweep > / < v' e v' z >. 

v z velocity fluctuations [3,4]. Following the def- 
initions given in [5] in a fixed frame, a strong 
sweep is associated with —v' 9 v' z > f3v' e rms v' z rms 
and v' z < (quadrant Q4) and a strong ejec- 
tion With -V' 6 V' Z > /H.rmX.mu and v 'z > 

(quadrant Q2). In the first quadrant Ql, v' e > 
and v' z > 0, while in the third quadrant Q3, 
v ' e < and v z < 0. The quadrant analysis is 
applied at z + = 17 corresponding to the lo- 
cation of the maximum value of the turbulent 
shear stress. We have also considered the value 
P — 2 to determine strong events, as used in 
[3,4]. We display in figures 3a and 3b the varia- 
tions with Ar + = r ± r + (r + the wall coordinate 
in the radial direction) of the conditionally aver- 



aged Reynolds shear stress normalized by the un- 
conditionally ensemble averaged Reynolds shear 
stress < v' e v' z > near a strong ejection (fig. 3a) 
and a strong sweep (fig. 3b). The contributions 
of each quadrant are also presented. The figures 
3a and 3b clearly show that the ejection (Q2) 
and sweep (Q4) quadrants contribute much more 
to the Reynolds shear stress production than the 
two other quadrants. On the other hand, it seems 
that the weakness of the turbulence in the present 
simulation accentuates the features observed in 
previous works. The results obtained support the 
conclusions of Littell and Eaton [3] and Lygren 
and Andersson [5], in contrast with the findings 
of Kang et al. [4]: the asymmetries observed 
are dominated by Reynolds stress-producing co- 
herent structures (sweep and ejection). Lygren 
and Andersson [5] concluded that clockwise vor- 
tices contribute much more to the Reynolds shear 
stress than counter-clockwise vortices. The same 
behavior applies in the presence of a sweep event. 
In this case, the levels of the surrounding ejec- 
tions approach the strong sweep level and are 
even slightly beyond the fixed criterion condition 
(3 = 2, as seen in figure 3b, while the levels of 
sweeps around a strong ejection are less impor- 
tant (fig. 3a), in agreement with the results of [5]. 
Case 1 vortices, having induced near-wall veloc- 
ity in the direction of the crossflow, are found 
to be the major source of generation of special 
strong events. 

Conclusion 



DNS calculations have been used to describe the 
turbulent boundary layer along the stator side 
in a large enclosed rotor-stator cavity. For the 
rotational Reynolds number under consideration 
Re = 9.5 x 10 4 , the Bodewadt layer is 3D tur- 
bulent, whereas the Ekman layer on the rotor is 
still laminar. The transition to turbulence is asso- 
ciated with the onset of localized turbulent spots 
spiraling up along the stator side. The turbulence 



kinetic energy budgets have revealed that produc- 
tion is the major contribution with a maximum 
obtained at z + ~ 12 independently of the ra- 
dial location, confirming the self- similar behavior 
of the Bodewadt layer. The results of the quad- 
rant analysis support the conclusions proposed 
by [3,5]. The asymmetries observed by these au- 
thors have been clearly detected and the analysis 
of conditionally averaged streamwise and wall- 
normal velocity components confirms that these 
asymmetries mainly arise from the contributions 
of quadrants Q2 and Q4, responsible for the gen- 
eration of ejection and sweep events. Therefore, 
Case 1 vortices are found to be the major source 
of generation of special strong events. 
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